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Abstract 

We compute the differential momentum correlation function for hadrons produced 
in the decay of S-mesons. This measure of hadronic event shapes tests the free 6-quark 
decay picture for nonleptonic and semileptonic decays of S- mesons. Our results can also 
be applied to Bs and decay. 
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Free heavy quark decay is expected to describe accurately the inclusive weak decays 
of hadrons containing a heavy quark. There are two physical reasons for this. First, while 
the individual exclusive modes depend on the details of hadronization, inclusive quantities 
which sum over all possible decay modes are expected to be insensitive to these details, 
since for heavy quarks there are many states available to the decay products. Second, for 
large heavy quark mass mq, the heavy quark in the hadron is almost on-shell and moving 
with the same four- velocity as the hadron. It is therefore reasonable to treat its decay as 
that of a free quark. 

For semileptonic (and radiative) decays, this duality between hadrons and free quarks 
and gluons has been shown to follow from QCD as the leading term in an expansion in 
l/mg Ij^. At leading order the differential decay rate for inclusive B — > XeVe decay, 
dV / dq^dEf, {q = + p,y), suitably averaged over E^,, is equal to that for free 6-quark 
decay. Non-perturbative corrections, parameterising the effects of the strong interactions 
of the heavy 6-quark with the light degrees of freedom in the S-meson, first arise at 
order l/m^. These have recently been calculated 0-0- Duality in this instance arises 
because the current-current correlator for the weak hadronic current is analytic everywhere 
in the complex energy plane, except at points on the real axis corresponding to physical 
intermediate states. This allows the phase space integral over the energy of the final 
hadronic state to be deformed to lie far from the the singularities in the decay amplitude. 
Along the deformed contour perturbative QCD provides a valid description of the decay 
products free of infrared singularities. The free quark decay picture emerges at leading 
order in 1/mf,. 

For nonleptonic decays the above argument does not apply, since there is no longer 
a kinematic variable that can be used to define a path of integration to be deformed. 
Nevertheless, the free quark decay picture is physically reasonable and widely used (e.g., 
to extract |Vc6| from the S- meson lifetime). Even if the free quark decay picture holds 
when nib is very large for the physical value of the 6-quark mass there is only about 
mb — rric — 3 GeV of energy (beyond the charm quark mass) available to the final state 
particles in weak B decay. Furthermore, this energy is shared amongst three quarks in 
b — > cud nonleptonic decay so that the energy that each quark carries is not very large 
compared with the QCD scale. It is important to test the free quark decay picture in as 
many ways as possible. One test is provided by the predicted equality of the S*^, B~ , Bs 
and Ab lifetimes and another by the predicted semileptonic branching ratios. In this 
letter, we point out that event shapes, characterized by the momentum angular correlation 
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distribution, may also be used to test its validity. We calculate this correlation function 
for both nonleptonic and semileptonic b decays assuming that the free quark decay picture 
is valid. 

Consider the weak decay of a hadron containing a 6-quark. We define the follow- 
ing quantity, in analogy with the differential energy-energy correlation function used to 
describe hadronic events at colliders [Q]: 



r is the total nonleptonic or semileptonic decay width, rrn, is the mass of the initial heavy 
6-quark, and the double sum runs over all pairs of strongly interacting particles (i, j), i j, 
in the final state separated by an angle 9ij, treating and (j, as distinct pairs. As 
we have defined it, the correlation function (|l]) is insensitive to the emission of soft and 
collinear gluons by the partons in the final state, which is required for it to be free of 
infrared divergences. Consequently dT,/d cos 9 at the hadron level should be equal (up 
to corrections of order 1/mf,) to its value calculated from free 6-quark decay. For b —>■ c 
transitions the final state contains a charmed hadron which decays weakly. It is important 
in the evaluation of (1) that this charmed hadron be treated as a single particle. One 
should not evaluate the sums in (1) using the charmed hadron's decay products. This 
complication makes dT, / d cos 9 more difficult to measure. 

For nonleptonic b decay, at leading order in a^, the tree level b qud matrix element 




(1) 
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of the weak Hamiltonian gives 



3Gj,ml [\pq\y/{yq - \pq\ cos 6 gd) (l - Vq - y a) 



dyq dya d cos 0qd 2ti^ 



1 - 2yg + ml 

2(1 -yq + \pq\ COs9qd) 



d^r 



3G%ml ( \Pq\yiP{yq - ml- yjtyq + yu\Pq\ COS Oqu 



dy q dy ud COS 9 qu 2tx^ 



X 5 h/ir- 



l-2yq + ml 



(2) 



d^V 



dy u dy dd cos 9 ud 



2(1 -yq+ \Pq\ COs9qu) 

3G|m| f yTc^yd^O- -yu + j/u cos 9ud) 
1 - yu + yu cos 9ud 
l-2yu- ml 
2(1 - yu + yu cos 9ud 
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where yq = Eq/mt, yu = E^/mb, yd = Ed/mb, Pq = Pq/mb and rhq = mq/mb- For the 
case, q = u, the energy integral in (|ip may be performed analytically, and the correlation 
function is found to be 



dT. 



dcos 9 



mq=0 



{l-cos9y 



(1 - COS 9) (811 + 1163 cos 9 + 401 cos^ 9 + 25 cos^ 9) 



+ 6(1 + cos ^) (195 + 171 cos ^ + 33 cos^ 9 + cos^ 9) log 
+ il5(l_cos^). 



1 + cos 6* 
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When q = c the mass of the charm quark cannot be neglected and, the phase space 
integration must be done numerically. In fig. |1| the momentum correlation functions for 
q = u (solid curve) and q = c (dashed curve) are shown for = 4.8 GeV and mc = 
1.6 GeV. In both cases, the correlation function is strongly peaked at ^ = tt, along with 
the delta function at 6* = 0. In addition to corrections suppressed by powers of l/mb, there 
are computable O (as(mb)) corrections to the results in fig. ^ 

For comparison with experiment, it is also useful to define the Fox- Wolfram moments 
H{L) [| of the distribution (|lD: 



H{L) 



dcos 9 



dT. 
dcos 9 



Pl{cos9) 
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where Pl{x) are Legendre polynomials. The first few Fox- Wolfram moments H{L) — H{0) 
(we have subtracted H{0) to remove the effect of the delta function at 6* = 0) are presented 
in Table 1. H{1) is equal to zero because of momentum conservation. Conservation of 
energy implies, in the massless q = u case, that -ff(O) = 1. 

We now consider inclusive semileptonic B decays. At leading order in as there is only 
one quark in the final state, so the correlation function is proportional to a delta function 
at ^ = 0: 

(iE(°) 2 1 - 9m^ + 45m^ - 45m^ + 9m^° - iv}^ + 120m^ log rhg 

dcos6 15 1 — Srhg + 8m^ — — 24mg log rfig (^5^ 

X 5(1 — cos 9). 

The leading perturbative contribution to the semileptonic correlation function away from 
cos^ = 1 arises from gluon bremsstrahlung at 0{as), and may be calculated using the 
results of |]T^. The rate for the process h — > qgeVe is 

1 1 f d?pc cPk d^pe d^p^ 4 1 



f d^Pc d'^k d'^ped^p^ . ,A,..,2 



2mb (27r) 

where the spin averaged matrix element for is 

\\M\^ = GWsTr \i>el^i>.lA\ - 75)] ■ 

rV Pc-k p,-k j^'^P p^-k p,.k ^')_ ' 

and k is the gluon four-momentum. The integration of the lepton phase space gives 

/ [l^e7M^7/3(l - 75)] = y [P,Pp - 9,,P'] , (8) 

where P = Pf, + p^- Once again, the integration can be performed analytically when the 
final quark is massless, and the 0{as) contribution to the correlation function is found to 
be 



dcos9 



4a^ /(83-5cos^)(l + cos^)3 l + cos^ 

log 



m -n Stt V 3(1- cos 2 

6875 + 17373 cos 9 + 12168 cos^ 9 + 1468 cos^ 9 - 507 cos^ ^ + 63 cos^ 9 ^ 
^ 360(1 - cos 

where we have neglected the 0{as) contribution to the delta function. The contribution of 
dT,^^^ / d cos 9 , in eq. (9), to H{0) diverges. This divergence is cancelled by contribution of 
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the order part of the delta function (that we neglected) to H{0). For a massive c quark 
in the final state, we have performed the phase space integration numerically. The results 
are plotted for q = u and q = c in fig. ||. Unfortunately, for both q = u and q = c there 
is only a small contribution to the correlation function away from the forward direction, 
making these distributions extremely difficult to measure. 

The first few Fox- Wolfram moments H[L) — H{0) for semileptonic b decays are pre- 
sented in Table 2. Note that these moments are much larger for b u transitions 
than b ^ c transitions. The large charm quark mass suppresses the amplitude for gluon 
bremsstrahlung. 

Clearly the results given in this paper also hold for A;, and Bs decay. For the q = u 
case, where the final quark is massless, our results can be applied to charm decay. However, 
the difference in lifetimes between the and is a clear sign that the free quark decay 
picture breaks down badly for nonleptonic charm decays, and so we do not expect our 
results to be useful in that case. 

Finally, we note that if the predictions for b decay event shapes that we have made 
disagree with experiment, it will not provide conclusive evidence that the free quark decay 
picture fails for the total nonleptonic decay width. It is possible that corrections suppressed 
by powers of l/mt are larger for event shapes than for the nonleptonic width. 

We thank H.D. Politzer and L. Wolfenstein for useful discussions. This research was 
supported in part by the Department of Energy under contract DE-FG02-91ER40682 and 
DE-FG03-92-ER40701. 
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Figure Captions 



Fig. 1. The momentum correlation function dTi/d cos 9 for nonleptonic decays computed 
in the parton model to lowest order in ctg. The solid curve corresponds to the 
case when all three quarks in the final state are massless. The dashed curve 
corresponds to the case for a h quark decaying to a charmed final state. 

Fig. 2. The momentum correlation function for (a) h ugeVf, and (b) b cgeVg as 
a function of cos^, omitting the normalization factor 4Q:s/37r. The dominant 
contribution to the delta function at ^ = comes from the tree-level graph; away 
from ^ = the distribution arises from the parton-level decay b — > qgeVf.. 
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H{L) 


- H{0) 


L 


ruq = 


rUq = 1.6 GeV 


1 


-1.0 


-0.72 


2 


-0.44 


-0.31 


3 


-0.74 


-0.54 


4 


-0.57 


-0.41 


5 


-0.67 


-0.49 


6 


-0.61 


-0.44 



Table 1: Fox- Wolfram Moments for Nonleptonic B Decay 





37r/4as(m6) 


x{H{L)-H{0)) 


L 


rUq = 


rUq = 1.6 GeV 


1 


-0.15 


-0.028 


2 


-0.26 


-0.037 


3 


-0.32 


-0.035 


4 


-0.36 


-0.035 


5 


-0.40 


-0.034 


6 


-0.43 


-0.034 



Table 2: Fox- Wolfram Moments for Semileptonic B Decay 
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